function f = parabola(N)
x = linspace(0,1,N);
f = zeros(N,1);
for i=1:N
    f(i) = 4*x(i)*(1-x(i));
end
end
